Slip distribution, coseismic deformation and Coulomb stress change for the 12 May 2008Wenchuan (China, Mw7.9) earthquake

The May 12, 2008 Wenchuan earthquake (Mw7.9) took place at the transition between the mountainous chain of Shan and the basin of Sichuan along the Longmen Shan Fault zone (31.1oN, 103.3oE; USGS). With a magnitude of 7.9 and a depth of ∼19 km the earthquake produced a 300-km-long fault rupture. It was the largest earthquake recorded in the region during the last centuries. It claimed more than 69,000 lives, induced widespread destruction over the region and raised concern about seismic hazard and source characterization for the Sichuan province. In the frame of our study, we selected 40 broadband waveforms (IRIS Consortium, USA) with good quality and satisfactory azimuthal coverage. Body waveforms were prepared for inversion using Kikuchi and Kanamori’s method [1] to obtain the spatiotem- poral slip distribution of a finite rupture model (length=300 km, strike=229o, dip=33o, width=60 km). The slip distribution model obtained was used to determine the coseismic deformation and the stress change distribution using the Coulomb 3.0 software [2]. Our coseismic deformation results was compared with data from GPS stations located near the fault rupture. Results show that directions of coseismic deformations are consistent with GPS observations close to the fault. Finally, we compare aftershock hypocenters that occurred during one month after the main shock with the Coulomb stress changes caused by this shock in the region. We observed that most aftershocks are located along the main fault plane without any noticeable clustering in the areas of increased stress. Our results suggest the rupture of the 2008 Wenchuan earthquake was essentially unilateral, from SW to NE (N49E), covering a 260km length and with duration about 105 sec. The strongest moment release occurred about 85km from the hypocenter, ∼30sec after the start of the rupture. Motions are dominated by thrust mechanism, but the superficial section of the second half of the rupture also shows a significant strike-slip component. [1]- Kikuchi, M., and Kanamori, H., 1982, Inversion of complex body waves: Bull. Seismol. Soc. Am., v. 72, p. 491-506. [2] -King, G. C. P., Stein, R. S. y Lin, J, 1994, Static stress changes and the triggering of earthquakes. Bull. Seismol. Soc. Am. 84,935-953

Influence of model parameters on synthesized high-frequency strong-motion waveforms

Waveform modeling is an important and helpful instrument of modern seismology that may provide valuable information. However, synthesizing seismograms requires to define many parameters, which differently affect the final result. Such parameters may be: the design of the grid, the structure model, the source time functions, the source mechanism, the rupture velocity. Variations in parameters may produce significantly different seismograms. We synthesize seismograms from a hypothetical earthquake and numerically estimate the influence of some of the used parameters. Firstly, we present the results for high-frequency near-fault waveforms obtained from defined model by changing tested parameters. Secondly, we present the results of a quantitative comparison of contributions from certain parameters on synthetic waveforms by using misfit criteria. For the synthesis of waveforms we used 2D/3D elastic finite-difference wave propagation code E3D [1] based on the elastodynamic formulation of the wave equation on a staggered grid. This code gave us the opportunity to perform all needed manipulations using a computer cluster. To assess the obtained results, we use misfit criteria [2] where seismograms are compared in time-frequency and phase by applying a continuous wavelet transform to the seismic signal. [1] - Larsen, S. and C.A. Schultz (1995). ELAS3D: 2D/3D elastic finite-difference wave propagation code, Technical Report No. UCRL-MA-121792, 19 pp. [2] - Kristekova, M., Kristek, J., Moczo, P., Day, S.M., 2006. Misfit criteria for quantitative comparison of seismograms. Bul. of Seis. Soc. of Am. 96(5), 1836–1850

Seismic source in the Iberian-African plate boundary

The plate boundary between Iberia and Africa has been studied using data on seismicity and focal mechanisms. The region has been divided into three areas: A; the Gulf of Cadiz; B, the Betics, Alboran Sea and northern Morocco; and C, Algeria. Seismicity shows a complex behavior, large shallow earthquakes (h < 30 km) occur in areas A and C and moderate shocks in area B; intermediate-depth activity (30 < h < 150 km) is located in the depth earthquakes (h » 650 km) are located to the south of Granada. Moment rate, slip velocity and b values have been estimated for shallow shocks, and show similar characteristics for the Gulf of Cadiz and Algeria, and quite different ones for the central region. Focal mechanisms of 80 selected shallow earthquakes (8 ‡ mb ‡ 4) show thrust faulting in the Gulf of Cadiz and Algeria with horizontal NNW-SSE compression, and normal faulting in the Alboran Sea with E-W extension. Focal mechanisms of 26 intermediate-depth earthquakes in the Alboran Sea display vertical motions, with a predominant plane trending E-W. Solutions for very deep shocks correspond to vertical dip-slip along N-S trends. Frohlich diagrams and seismic moment tensors show different behavior in the Gulf of Cadiz, Betic-Alboran Sea and northern Morocco, and northern Algeria for shallow events. The stress pattern of intermediate-depth and very deep earthquakes has different directions: vertical extension in the NW-SE direction for intermediate depth earthquakes, and tension and pressure axes dipping about 45 ° for very deep earthquakes. Regional stress pattern may result from the collision between the African plate and Iberia, with extension and subduction of lithospheric material in the Alboran Sea at intermediate depth. The very deep seismicity may be correlated with older subduction processes

Tomographic three-dimensional seismic velocity structure of the SW Ibero-Maghrebian region

The present tomographic study focuses on SW Ibero-Maghrebian region. To locate the seismic events and find the local velocity structure of epicentral area, the P and S arrivals at 42 stations located at north of Morocco, south of Portugal and Spain are used. The arrival times data used, in this study, were obtained by the “Instituto de Meteorologia” (IM, Lisbon, Portugal), the National Institute of Geophysics (CNRST, Rabat, Morocco) and the “Instituto Geografico Nacional” (IGN, Madrid, Spain) (between 12/1988 and 30/2008). The preliminary estimate of origin times and hypocentral coordinates are determined by the hypocenter 3.2 program. In this study we use a linearized inversion procedure comprising two steps: 1) finding the minimal 1-D model and simultaneous reloca- tion of hypocenters and 2) determination of local velocity structure assuming a continuous velocity field. The earth structure is represented in three dimensions by velocity at discrete points, and velocity at any intervening point is determined by linear interpolation among the surrounding eight grid points. The resolutions tests results indicate that the calculated images give near true structure for the studied region at 15, 30, 45 and 60 km depth. At 5km depth it gives near true structure in the continental region of Portugal, Spain, and Morocco. This study shows that the total crustal thickness varies from 30 to 35 km and contains low-velocity anomalies. A prominent low velocity anomaly that shows a maximum decrease in P-wave velocity of approximately 6 per cent in the Gibraltar region is observed extending down to a depth of approximately 30 km. This low velocity demarcates a small bloc located between Iberia and Nubia plates. The resulting tomographic image has a prominent high velocity anomaly that shows a maximum increase in P-wave velocity of approximately 6 per cent between 45 to 60 km depth beneath South of Portugal and the Golf of Cadiz. High-velocity anomalies could be associated with the location of deep active faults in the uplift and upper crust of South of Portugal. In the Golf of Cadiz, these anomalies could be associated with the seismogenic zone and probably more at the south with the Iberia-Nubia plate boundary

Design and Evaluation of a High Throughput Seismic Sensor Network - Tools for Planning, Deployment and Assessment

The rapid technological evolution in sensors, sensor platforms and networking is enabling the deployment of large sensor networks for "live" monitoring of seismic activity with high spatial resolution. In this regard, this paper describes our work in developing an online "High Throughput Seismic Sensor Network". We present the architecture and implementation comprising seismic sensors and servers (running data collection services) connected through internet-enabled technologies. We validate and assess the system, as well as identify bottlenecks, by means of experimentation. Based on the collected empirical data, we were able to identify methods and tools to support effective planning and implementation of sensor networks based on two main indicators: Sensor Network Transmission Rate (SNTR), which provides the overall network sensor data transmission throughput and thus an indication of the required network capacity; and CPU Sensor Network Performance Index (CSNPI), which provides an indication of a server capability to handle network sensor data. As we progress in our work to field deploy seismic sensor networks, we will continue to use these tools to plan and deploy future sensor networks, as well as assess improvements and modifications along the way

PLASMA – a high-performing and open platform for the integration of heterogeneous sensor networks

The use of sensors to capture disparate types of information from the environment has been increasing and they cover a wide range of applications, such as climate monitoring (e.g., seismic activity and climate), water quality monitoring, area surveillance, intelligent buildings, energy management, automotive industry and scientific purposes. Additionally, the Information Age has provided rapid and ubiquitous access to information produced by heterogeneous sources. Thus, the development of sensor networks has emerged as a way to exploit the Information Age capabilities into sensor applications. A major outcome of this combination is the capability to remotely receive and process sensor data (covering large areas) in real-time therefore allowing the development of new studies and algorithms aiming at anticipating and predicting (with high- reliability) events, such as storms (already state-of-the-art) and earthquakes (not yet possible). The developments in the scientific and industry communities were proficient in creating a large number of sensor networks that, nonetheless, (i) were designed to fit specific needs, (ii) were usually deployed in closed networks (not accessible to external parties), and/or (iii) do not interoperate with each-other, thus not leveraging possible synergistic effects of combining multiple sensor networks. To solve these issues, we initiated the project PLASMA

DIRDOP: a directivity approach to determining the seismic rupture velocity vector

Directivity effects are a characteristic of seismic source finiteness and are a consequence of the rupture spread in preferential directions. These effects are manifested through seismic spectral deviations as a function of the observation location. The directivity by Doppler effect method permits estimation of the directions and rupture velocities, beginning from the duration of common pulses, which are identified in waveforms or relative source time functions. The general model of directivity that supports the method presented here is a Doppler analysis based on a kinematic source model of rupture (Haskell, Bull Seismol Soc Am 54:1811–1841, 1964) and a structural medium with spherical symmetry. To evaluate its performance, we subjected the method to a series of tests with synthetic data obtained from ten typical seismic ruptures. The experimental conditions studied correspond with scenarios of simple and complex, unilaterally and bilaterally extended ruptures with different mechanisms and datasets with different levels of azimuthal coverage. The obtained results generally agree with the expected values. We also present four real case studies, applying the method to the following earthquakes: Arequipa, Peru (Mw = 8.4, June 23, 2001); Denali, AK, USA (Mw = 7.8; November 3, 2002); Zemmouri– Boumerdes, Algeria (Mw = 6.8, May 21, 2003); and Sumatra, Indonesia (Mw = 9.3, December 26, 2004). The results obtained from the dataset of the four earthquakes agreed, in general, with the values presented by other authors using different methods and data

Strong Ground Motion Simulations and Assessment of Influence of Model Parameters on the Waveforms

Modeling near-field ground motion is an important and useful tool of modern seismology. In our work we use a finite difference algorithm to compute near-field ground motions from a real moderate event with pre-existing slip distribution model. Lately, synthetic seismograms are quantitatively compared with observed waveforms from near-field seismic stations in order to justify created model. Furthermore, we independently changed several source parameters (rupture velocity, source dimension and geometry), and structure (velocity model) in order to evaluate their influence on the waveforms. For the comparison of seismograms we applied quantitative misfit criteria based on wavelet transform

Mass – Radius Relationship in Extrasolar Planets

The increasing number of Extrasolar planets observed in the last years makes important to define, as soon as possible, a mass – radius relationship, and so, we adjusted an planetary constitution independent experimental equation. Using the latest database of Extrasolar planets, a bi-logarithmic graphic was plotted that represents the mass - radius relationship where we adjusted a polynomial equation, which better suited the sample of Extrasolar planets at current time